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Abstract 

We study the stress distribution profiles and the height 
and thickness fluctuations of lipid membranes in the 
gel Lpi state by Monte Carlo simulations of a generic 
coarse-grained model for lipid membranes, which re- 
produces many known properties of dipalmitoylphos- 
phatidyncholine (DPPC) bilayers. The results are re- 
lated to the corresponding properties of fluid mem- 
branes, and to theoretical predictions for crystalline and 
hexatic membranes. One striking observation is that 
the spontaneous curvature of the monolayers changes 
sign from the fluid to the gel phase. In the gel-phase, 
the long-wavelength height fluctuations are suppressed, 
and the fluctuation spectrum is highly anisotropic. In 
the direction of tilt, it carries the signature of soft modes 
that are compatible with the wavelength of the ripple 
phase Ppi, which emerges in the transition region be- 
tween the fluid and the gel state. In the direction per- 
pendicular to the tilt, the thickness fluctuations are al- 
most entirely suppressed, and the height fluctuations 
seem to be dominated by an interfacial energy, i.e., by 
out-of-layer fluctuations, up to length scales of tens of 
nanometers. 

1 Introduction 

Lipid molecules in water self-assemble into a variety 
of nanostructures. Among the most prominent are bi- 
layers, where the lipids are arranged such that their hy- 
drophobic hydrocarbon tails are shielded from the water 
by the hydrophilic head groups. Lipid bilayers form the 
matrix for biological membranes, and they are ubiqui- 
tous in all living matter . 



Model systems of pure (one-component) lipid mem- 
branes in water have been studied extensively for many 
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decades — — . Depending on the temperature and on the 
structure of the lipids, one finds a variety of phases: 
A high-temperature "fluid" phase (L a ) where the lipid 
tails are highly disordered, a low-temperature "gel" 
phase with either straight (Lp) or collectively tilted 
(Lpi) tails, and in the latter case an intermediate modu- 
lated "ripple" phase (Pp/). At still lower temperatures, 
other even more ordered phases can be identified. Na- 
ture usually maintains biomembranes in the fluid state. 
Nevertheless, it is remarkable that the transition be- 
tween the fluid and the gel state, the so-called main 
transition, occurs at temperatures well above room tem- 
perature for some of the most common lipid molecules, 
(e.g., ~ 41°C for the phospholipid dipalmitoylphos- 

phatidyncholine (DPPC)) ~. 

The gel state of membranes differs from the fluid 
state in several respect. As mentioned above, the hy- 
drocarbon chains are much more ordered i.e., they 
have fewer gauche bonds. As a result, the membrane 
is thicker and the nematic order parameter inside the 
membrane is higher. Furthermore, the molecules also 
have a higher degree of (two-dimensional) positional 
order within the membrane: Gel membranes are pre- 
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sumably hexatic or (quasi)crystalline . Both factors, 
the higher chain order and the higher positional order, 
significantly affect the elastic properties of membranes: 
They are much stiffer in the gel state than in the fluid 
state. 

The simplest Ansatz for the curvature elastic free en- 
ergy density per unit area of a tensionless fluid mem- 
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brane is given by 

,^/A = 2K{H-c ) 2 + kK, (1) 

where 2H is the total curvature of the membrane, i.e., 
the sum of the inverse curvature radii, and K is the 
Gaussian curvature, i.e., their product. The elastic pa- 
rameters K, k, and cq are the bending rigidity, the Gaus- 
sian rigidity, and the spontaneous curvature of the mem- 
brane.- For symmetric bilayers, cq is zero. Further- 
more, the Gauss-Bonnet theorem states that the inte- 
gral over K on closed membranes contributes a constant 
which only depends on the topology, hence it can be 
omitted if the latter is fixed. The only remaining pa- 
rameter is the bending rigidity K. 

In hexatic or crystalline membranes, the situation is 
more complex due to the fact that lipids have fixed 
neighbors and membranes can sustain elastic shear 
stress. In a seminal paper of 1987 , Nelson and Peliti 
have pointed out that in this case, fluctuations stiffen the 
membranes, and the effective bending rigidity appears 
to be increased on larger length scales. Out-of-plane 
thermal undulations with different wavevectors q are 
accompanied by in-plane shear deformations and thus 
penalized. This leads to long-range phonon-mediated 
interactions between capillary waves, which renormal- 
ize the free energy For crystalline membranes, Nel- 
son and Peliti predicted that the effective bending rigid- 
ity scales as Kft(q) l/q as a function of the wavevec- 
tor. Later, this estimate was refined by Le Doussal and 

Radzihovsky , who obtained fCff(g) « q~ S21 within 
a self-consistent screening approximation. In hexatic 
membranes, the positional order is destroyed due to the 
presence of free dislocations, but they still have long 
range 'bond-orientational' order, i.e., the vectors con- 
necting nearest neighbor lipids have well-defined aver- 
age orientations. According to Nelson and Peliti, the 
bending stiffness then still increases for large length 
scales, but only logarithmically, KR(q) ~ yj — ln^^y), 
where §y is a translational correlation length. Nelson 
and Peliti also argued that physical lipid membranes 
should be hexatic rather than crystalline, because the 
ability to buckle reduces the elastic energy of free dis- 
locations. 

The situation is further complicated in the Lp> phase 
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due to the collective tilt of the molecules. Park "" has 
studied fluctuating tilted hexatic membranes within a 
model where two order parameter fields, one describ- 
ing tilt order and one bond-orientational order, evolve 
on a fluctuating fluid surface ([]]). Local shear modes, 



which play a major role in the bending renormalization 
mechanism of Nelson and Peliti, are not included. In 
this theory, the bending rigidity in the so-called 'strong 
coupling limit' renormalizes to a fixed ^-independent 
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value, both for untilted and tilted membranes. 

It should be noted that fluctuations also renormal- 

ize the bending rigidity in fluid membranes, albeit in 

the opposite way than in hexatic/crystalline membranes. 

On large length scales the membranes become softer, 

and Kr effectively vanishes in the infinite wavelength 
n 

limit . However, this softening becomes only relevant 
on length scales where the membranes exhibit large de- 
viations from planar. In practical simulation studies 
of fluctuating self-assembled planar fluid membranes, 
a well-defined bending rigidity that can be associated 
with the unrenormalized (local) bending rigidity can 
usually be extracted from the capillary fluctuation spec- 
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trum - — . 

The purpose of the present work was to study the 
elastic properties of membranes in the fluid and the gel 
state, and in particular, to look for signatures of the 
liquid/gel transition in the membrane fluctuation spec- 
trum. To this end we use a coarse-grained approach to 
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membrane modeling . We have carried out exten- 
sive simulations of a very simple generic lipid model 
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, which has been shown to self-assemble into bi- 
layers with reasonably realistic elastic parameters in 
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the fluid phase and to reproduce the most impor- 
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tant membrane phase transitions of DPPC , i.e., the 
transition between the L a -phase, the Ppi phase, and the 
Lpi -phase. In previous work, we have characterized 
the fluid phase and analyzed the resulting membrane- 
protein and membrane-mediated protein-protein inter- 
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actions — . Now, we will consider the Lpi phase, fo- 
cussing on the stress profiles in the membranes, which 
can be related to the elastic parameters, and on the cap- 
illary fluctuations. 

The paper is organized as follows: In the next section, 
we introduce the model and the simulation method. The 
results are presented in section[3] Finally we summarize 
and conclude in section [4] 

2 Model and Methods 

In our model, each lipid is represented by a chain of 
seven beads with one head bead of diameter (7/, fol- 
lowed by six tail beads of diameter a,. Non-bonded 
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beads interact, via a truncated and lifted Lennard-Jones 
potential: 



V M (r) 



with 



V u {r/a)-V u {r c /a) if r<r c 
otherwise 



V u 0) = e h 
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2x 
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(2) 



(3) 



where a is the mean diameter of the two interacting 
beads, (7;; = (<7; + (7;)/2 = /i or f). Head-head 
and head-tail interactions are purely repulsive (r c = (7), 
while tail-tail interactions also have an attractive con- 
tribution (r c = 2(7). Bonded beads are connected by 
FENE (Finitely Extensible Nonlinear Elastic) springs 
with the spring potential 

V FENE (r) = -ie FENE (Ar max ) 2 log ^1 - 

' (4) 

where ro is the equilibrium distance, Ar max the maximal 
deviation, and £ FENE the FENE spring constant. In addi- 
tion, we introduce a bond-angle potential 




V B a(0) = £ba(1-cos(0)). 



(5) 



The aqueous environment of the membrane is modeled 
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with "phantom" solvent beads , which interact with 
lipids like head beads (a s = (7/,), but have no interac- 
tions with each other. 
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The model parameters are " o), = 1.1(7,, ro — 
0.7(7,, Ar m „ = 0.2(7,, = 100e/(7, 2 , and e BA = 4.7e, 
and the pressure was chosen P — 2.e/a, 3 . At these pa- 
rameters, the lipids spontaneously self-assemble into 
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stable bilayers , and the bilayer undergoes a main 
transition to a tilted gel phase Lpi via a ripple phase 
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Ppi at the temperature ksT = 1 .2e 

The system was studied using Monte Carlo 
simulations at constant pressure and temperature 
with periodic boundary conditions in a simulation 
box of variable size and shape: The simulation 
box is a parallelepiped spanned by the vectors 
(L*,0,0), (syxL^Ly^^SzxL^SzyLy,^), and all L, and 
sj are allowed to fluctuate. This ensures that the mem- 
branes have no interfacial tension. We have verified 
that the total pressure tensor had the form £? a $ = Pd a p 
with the applied pressure P. We have studied systems 
of up to 7200 lipids, and the run lengths were up to to 
8 million Monte Carlo steps, where one Monte Carlo 
step corresponds to one Monte Carlo move per bead, 



and moves that change the shape of the simulation box 
were attempted every 50th Monte Carlo step. The code 
was parallelized using a domain decomposition scheme 
described in Ref. 21 . The configurations were initially 
set up as an ordered bilayers with straight chains in the 
(xy)-plane. Typical equilibration times were 1 million 
Monte Carlo steps for the small systems, which we used 
to determine pressure profiles, and 4 million Monte 
Carlo steps for the large systems, which we used to ex- 
tract fluctuation spectra. After this time, the results for 
the fluctuation spectra did not change any more. In the 
fluid state, one can potentially reduce the equilibration 
and the correlation times by applying additional cluster 
moves that specifically speed up the long-wavelength 
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modes — . This was not done here, since these modes 
were not the most critical ones in the gel state (see be- 
low), and the length of the runs was mainly dictated by 
the need to accumulate good statistics over the whole 
range of wave vectors. It should be noted that in the 
gel state, lipids never flipflopped from one monolayer 
to the other, hence the number of lipids on both sides 
remained equal throughout the simulation. 

The model units can approximately be translated into 
standard SI units by comparing the bilayer thickness 
(~ 6(7, in the L a phase and ~ 7.7(7, in the Lpi phase) and 
the area per lipid (~ 1 .4a 2 in the L a phase and ~ 1 of in 
the Lpi phase) with the corresponding numbers for real 
lipid bilayers. The numbers in our model roughly repro- 
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duce those of DPPC bilayers if we identify (7, ~ 6A — . 
The energy scale can be estimated by matching the tem- 
peratures of the main transition: e ~ 0.36 • 10~ 20 J. 



3 Results 

The simulations were mainly carried out at kgT = 1 .0e, 
which is well in the Lpi phase. For comparison, we also 
show data for fluid membranes at IcbT = 1.3e. Fig. [T] 
shows slices through membrane configurations in the 
fluid and the gel state. The corresponding structure fac- 
tors (data not shown) feature distinct peaks in the gel 
state, corresponding to triangular positional order, and 
diffuse Debye-Scherrer rings in the fluid state, reflect- 
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ing lack of positional order in that phase — . At the 
available system sizes, it was not possible to determine 
whether the positional order in the gel phase is hexatic 
or crystalline. 

We begin with considering the local pressure profiles, 
which give information on the internal stress distribu- 
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Figure 1: Slices through configurations (side view) of a 
self-assembled lipid bilayer in the fluid phase (top) and 




Figure 2: Top left: Interfacial tension profile y(z) for 
lipid bilayers with N = 200 (solid line) and N = 3200 
(dashed line) in the gel phase. In large systems, the pro- 
files are slightly broadened due to the fluctuations of 
the membrane position. The subsequent profiles were 
evaluated for systems with N = 200 lipids. Bottom left: 
Corresponding density profiles of solvent (dotted line), 
head (dashed lines), and tail beads (solid lines). Right: 
Contributions to the surface tension profile. Top right: 
Contributions from nonbonded interactions (NB, thick 
solid line), bond angle potentials (BA, dotted line), and 
bond length potentials (BL, dashed line). Bottom right: 
Among the nonbonded interactions, contributions from 
head-head interactions (HH, solid line), head-tail inter- 
action (HT, dotted line), head-solvent interaction (HS, 
dashed line), and tail-tail interaction (TT, dot-dashed 
line). The contribution from tail-solvent interactions is 
negligible. 



tion in the membranes. The pressure tensor was deter- 
mined using a version of the virial theorem that contains 
only relative coordinates, 




where N is the number of particles, T the temperature, 
V the volume, ry = r, — rj the vector connecting the 
particles i and j, Fy :paiI the pairwise force acting on par- 
ticle i from particles /', and 0, the angle between the 
two bonds from bead j. To determine the local pres- 
sure profiles, the system was divided in bins of length 
8z = 0.125o> along the z-axis. The contributions of rel- 
ative vectors r,-; were parceled out on the bins according 
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to the Irving-Kirkwood convention — . 

Fig. |2] shows profiles of the local interfacial tension 
y(z), defined as the difference Pn(z) —Pt(z) of the nor- 
mal and the tangential pressure profile, 

7(2) = & a {z) - {Pa® + ^yy(z))A (7) 

for membranes in the gel state. We note that for our 
tensionless membranes, the integral J dz y(z) vanishes. 
The contributions of different forces to the full profile 
are shown in Fig. [2] (right). 

The main features of the profile can be summarized 
as follows: It features (i) maxima in the regions where 
head and tail beads meet, reflecting the interfacial en- 
ergy of the hydrophilic/hydrophobic interface, (ii) min- 
ima in the head regions, reflecting the tendency of the 
heads to spread out, and (iii) a minimum at the cen- 
ter of the bilayer, reflecting the 'negative interfacial 
tension' or adhesive coupling between the monolayers 
{i.e., a tendency of the monolayer-monolayer interface 
to spread out). These peaks were also found in the 
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stress profiles of fluid membranes in our " and sim- 
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ilar — coarse-grained models, as well as in atomistic 
and more reastically coarse-grained simulations of fluid 
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lipid bilayers . The additional weak maxima at the 
head/solvent interfaces reflect a head-solvent interfacial 
energy which is specific for our model and not always 
present. Overall, the stress profiles in the fluid and the 
gel phase are strikingly similar. They differ mainly in 
amplitude. Compared to the fluid state, the surface ten- 
sion profile is amplified by roughly a factor of ten in 
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the gel state, thus the membrane is under much higher 
internal tension. Furthermore, the relative depth of the 
minima in the head region increases, compared to the 
central minima, i.e., the spreading pressure in the head 
region gains importance. As a result, the monolayers 
develop a tendency to bend inwards, i.e. they acquire a 
spontaneous curvature. 

For fluid membranes, the latter can be estimated by 
the first moment of the stress profile,— kcq = —Zy with 

z r := f dz%*(z)z, (8) 
Jo 

where we have used / dz y(z) = 0, and fc/2 is the mono- 
layer bending rigidity. We emphasize that cq refers to 
the spontaneous curvature of a monolayer. The sponta- 
neous curvature of the bilayer is zero, as mentioned in 
the introduction. In the gel state, the interpretation of z.y 
is less obvious, since the Helfrich description (Q]) is no 
longer valid. As we shall show below, the very notion 
of a bending rigidity becomes questionable on nanome- 
ter length scales. Nevertheless, one can argue that iy 
still gives qualitative information on the propensity of 
monolayers to bend in one direction. 

Applying Eq. ® yields z r = -11.1 ± O.le/oj. In 
the fluid state, we had obtained z.y = +0.3 ± O.le/Cr 
in Ref.: 17 Thus fluid monolayers were found to have a 
weak tendency to bend outwards, in agreement with the 
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results from atomistic simulations ". Here, in the gel 
state, they have a tendency to bend inwards. We con- 
clude that the elastic properties of monolayers in the gel 
phase differ qualitatively from those in the fluid phase. 
The interplay of two spontaneous curvatures is presum- 
ably one of the factors that stabilize a modulated ripple 
phase in the transition region between the fluid and gel 
phase. 

Next we discuss the fluctuation spectra of the mem- 
branes. We have simulated lipid bilayers containing 
up to 7200 lipids, corresponding to patches with ar- 
eas of roughly 1300 nm 2 . The data were analyzed fol- 
lowing a procedure described in Refs. 15,17 . Basically, 
we determine the mean head positions zi(x,y), Z2(x,y) 
of the two monolayers on a grid (x,y), and derive the 
height and thickness profiles h(x,y) = (zi +zi)/2 and 
t(x,y) = (z\ —Z2}/2. The spectra are Fourier trans- 
formed according to 

U, 9 y = ^-Lf^y)e-^ x+w) , (9) 

and averages (\h q \ 2 ) and (|f 9 | 2 ) are evaluated in g-bins 
of size 0.1/ Of 
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Figure 3: Radially averaged Fourier spectrum of the 
height (full circles) and thickness (open squares) fluctu- 
ations of membranes in the gel state. Black solid lines 
show the best fit to the elastic theory by Brannigan and 
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Brown (see text for explanation). Grey solid line in- 
dicates the slopes corresponding to a g~ 4 -behavior and 
a g~ 3 -behavior. Inset (from Ref. 17 ) shows the corre- 
sponding data for the fluid membrane for comparison. 

Fig. |3] shows the radially averaged results. For com- 
parison, the corresponding data for the fluid phase 
(from Ref.—) are shown in the inset. In the fluid 
phase, the data could be fitted very nicely with an elas- 
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tic theory due to Brannigan and Brown , which de- 
scribes the bilayer in terms of two coupled elastic sheets 
(monolayers) with two types of independent fluctua- 
tions: "bending" fluctuations and "protrusion" fluctu- 
ations. The spectra are dominated by the protrusion 
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modes at high ^-vectors , and by the bending modes 
at low q- vectors. 

In the gel state, the qualitative features of the radially 
averaged spectra are at first sight very similar. For ex- 
ample, the height and the thickness spectra nearly fall 
on top of each other at high q, and the thickness spec- 
trum seems to feature a a soft peristaltic mode at q ~ 
0.5/(7 ( , corresponding to a wavelength around 10(7 f . 
Quantitatively, however, one notices distinct differences 
between the fluid and the gel state: First, the amplitudes 
of the fluctuations in the gel state are much smaller than 
in the fluid state. Second, the elastic theory fails at low 
g-vectors, even though it seems to describe the data rea- 
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Figure 4: Fourier spectrum of the height (full circles) 
and thickness (open squares) fluctuations of membranes 
in the gel state for wavevectors parallel (top) and per- 
pendicular (bottom) to the tilt direction. Thin vertical 
bar in the upper panel marks the position of a peak in 
both height and thickness fluctuations. Dashed lines in 
the lower panel show <7~ 2 -slope and g~ 4 -slope as indi- 
cated. 



sonably well at high g-vectors. Compared to the theo- 
retical prediction, the long-wavelength height fluctua- 
tions are suppressed. The low-g limit (\h(q)\ 2 ) ~ q~ 4 is 
never encountered. The data are more consistent with 
the asymptotic behavior predicted by Nelson and Peliti 

for crystalline membranes , (\h(q) \ 2 ) ~ 1/ K]t(q)q~ 4 ~ 
g~ 3 , but in a double-logarithmic plot, even the slope of 
q~ 3 is slightly too steep (Fig. 01. 

To explore this phenomenon in more detail, we have 
analyzed separately the fluctuations in tilt direction and 
perpendicular to the tilt direction. Since the tilt direc- 
tion slowly fluctuated throughout the simulations, this 
means that before averaging (\h(q)\ 2 ) and (|f(q)| 2 ), ev- 
ery configuration had to be rotated such that the tilt 
points along a given axis. The results are shown in Fig. 



|4] Both in the direction parallel and perpendicular to 
the tilt, the thickness and height fluctuation spectra dif- 
fer strongly from the radially averaged spectra, and are 
no longer compatible with the elastic theory. 

In the tilt direction, the thickness fluctuation am- 
plitudes are largely comparable to those of the radi- 
ally averaged spectrum, with the exception of the q- 
region around (q(7 t ) 2 ~ 0.15, where they exhibit a max- 
imum. In contrast, the height fluctuations are strongly 
suppressed, most notably for large g-vectors beyond 
(q<Jt) 2 ~ 0.7, but also for small g-vectors. Both the 
thickness and the height fluctuation spectra feature 
small but distinct peaks at qo t ~ 0.8, which clearly rise 
above the range of statistical fluctuations. These peaks 
correspond to undulations with wavelength A ~ 8(7,. 
The tendency of the gel membrane to exhibit structure 
on that length scale can already be sensed when looking 
at the bare snapshots of the gel phase (Fig.[T]i. The other 
maximum at (gcr,) 2 ~ 0.15 (or go 7 ,) ~ 0.4) in the thick- 
ness spectrum, which was mentioned in the beginning, 
corresponds to undulations with twice the wavelength, 
A ~ 16(7,. In the height spectrum, one can vaguely dis- 
cern a 'peak' in the same g-region, which however inter- 
feres with the general growth of fluctuation amplitudes 
for small g- vectors. 

Hence we conclude that both the thickness and the 
height fluctuation spectra feature soft modes with wave- 
length Ao ~ 8(7, and X\ = 2Xq. The wavelength Ai = 
16(7 happens to coincide with the period of the ripple 
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phase . This suggests that the fluctuations in tilt direc- 
tion already carry the signature of the modulated phase 
which intrudes between the liquid and the gel state. 

In the direction perpendicular to the tilt, the situation 
is different. The thickness fluctuations are almost en- 
tirely suppressed. The height fluctuations rise smoothly 
and are best compatible with a q~ 2 slope, suggesting 
that they are controlled by an interfacial tension rather 
than a bending rigidity. On the other hand, the evalu- 
ation of the pressure tensor clearly shows that the total 
interface tension, integrated over the membrane, van- 
ishes within the error (y=0±0.1e/<7, 2 ). 

An alternative potential explanation for the observed 
q~ 2 shape of the fluctuation spectrum can be derived 
based on the following two assumptions: (i) The mem- 
brane undulations in the direction perpendicular to the 
tilt (denoted x hereafter) are dominated by out-of-layer 
fluctuations, i.e., the lipids mainly move up and down, 
keeping their orientation, such that energy contributions 
due to bending or splay are negligible (see Fig- [5J. (ii) 
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Figure 5: Schematic sketch of different types of mem- 
brane height fluctuations. Left: Out-of-layer fluctua- 
tions, governed by a surface energy that depends on 
the local (projected) tilt angle 9. Right: Bending- 
dominated height fluctuations, governed by a bending 
energy. 

The effective interfacial tension in x-direction y x de- 
pends on the angle 9 between the projections onto the 
xz plane of the lipid orientation vector and the local 
interface normal. (We note that this second assump- 
tion implies that the x-direction can be defined inde- 
pendently of the tilt direction, e.g., in terms of the 
hexatic order in the membrane.) For small angles 9, 
the function % can then be expanded in powers of 9, 
Yx[8] w A9 2 /2. For constant lipid orientation, the an- 
gle 9 depends on the local membrane height position 
h(x,y) via 9 « dh/dx. Hence the fluctuations of h(x,y) 
in x-direction contribute to the total free energy with 
the term/ dxy x [9(x)} «A/2/ dx(dh/dx) 2 , which looks 
formally like an interfacial tension contribution with 
tension A. This explains how out-of-layer fluctuations 
can display a tension-like behavior. One would ex- 
pect the out-of-layer fluctuations to eventually give way 
to bending-dominated fluctuations at very long wave- 
lengths, but the latter are not observed in our simula- 
tions. 



4 Discussion and Summary 

To summarize, we have analyzed the elastic behavior 
and the fluctuations of lipid membranes in the gel state. 
They are distinctly different from those in the fluid state. 
Whereas the liquid state can be described nicely by an 
elastic theory for fluid coupled monolayers, this de- 
scription fails in the gel state. The failure can be at- 
tributed to the ordered structure of the gel membrane 
and to its anisotropy. Our simulation results disclose a 
complex fluctuation spectrum, which cannot be under- 
stood within any of the available continuum theories for 
membranes. 

The spectrum is anisotropic and features two soft 
modes that are very likely related to the ripple phase, 



which is known to emerge in the transition region be- 
tween the liquid and the gel phase. Given the fact that 
the latter has a periodic modulation with wavelength 
~ 10 nm, we must conclude that our systems, which 
have linear dimensions of around ~ 35 nm, are still 
much too small to give insight into the asymptotic elas- 
tic behavior of the membranes. Unfortunately, the sim- 
ulation of much larger systems (at least a factor of ten 
would be necessary) is unfeasible with our currently 
available computer resources. Our simulations can be 
used to analyze the elastic behavior of the membrane 
on the length scale of tens of nanometers. Based on 
our data, it is clear that the membrane in the gel state 
still has no single unique 'bending rigidity' on this lat- 
eral length scale, even though it exceeds the molecular 
scale (the monolayer thickness) by more than an order 
of magnitude. Rather, the elastic behavior results from a 
complex interplay of different contributions which can- 
not easily be told apart. Formulating an elastic theory 
for such membranes will be a challenge for future the- 
oretical work, which should presumably also contribute 
to a better understanding of the mechanisms driving the 
formation of the ripple phase. 

In the present work, we have discussed 'mon- 
odomain' gel membranes with a uniform tilt direction. 
This seems realistic, given that scanning optical micro- 
scope studies of supported DPPC-bilayers in the gel 
phase— have shown that single-tilt domains in such 
membranes extend over 1-2 micrometers. On larger 
scales, the membranes are effectively 'polycrystalline', 
and their fluctuations will be dominated by processes 
at the domain boundaries and/or defects in the tilt or- 
der. Preliminary studies, where we have enforced such 
defects, have shown that they amplify the fluctuations 
considerably. It will be interesting to study these effects 
in more detail in the future. 
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